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The operator-sum decomposition (OS) of a mapping from one density matrix to another has 
many applications in quantum information science. To this mapping there corresponds an affine 
map which provides a geometric description of the density matrix in terms of the polarization 
vector representation. This has been thoroughly explored for qubits since the components of the 
polarization vector are measurable quantities (corresponding to expectation values of Hermitian 
operators) and also because it enables the description of map domains geometrically. Here we 
extend the OS-amne map correspondence to qudits, briefly discuss general properties of the map, 
the form for particular important cases, and provide several explicit results for qutrit maps. We use 
the affine map and a singular- value-like decomposition, to find positivity constraints that provide a 
symmetry for small polarization vector magnitudes (states which are closer to the maximally mixed 
state) which is broken as the polarization vector increases in magnitude (a state becomes more 
pure). The dependence of this symmetry on the magnitude of the polarization vector implies the 
polar decomposition of the map can not be used as it can for the qubit case. However, it still leads 
us to a connection between positivity and purity for general d-state systems. 
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I. INTRODUCTION 

The study of maps of density operators has become an active field of research in quantum information theory. 
Such maps are used to describe open-system evolution, i.e., evolution of quantum systems which interact with their 
environment in a non-trivial way. Such maps are also useful for describing methods to alleviate such noise and, in 
some cases, as indicators of entanglement [l|, 0. These are especially important in the study of quantum information 
processing. 

General maps of density operators to density operators were first studied by Sudarshan, Mathews and Rau (SMR) 
Q. They provided a method for describing open-system evolution and referred to such maps as "dynamical maps." 
This name is particularly appropriate for describing noise since detailed properties of the environment causing the 
noise are often unknown. Often these maps are written in what is called an operator-sum representation (OSR). Kraus 
later considered completely positive maps [J. The assumption of complete positivity is sometimes useful but is not 
a necessary assumption for open system evolution. The physical implications of the complete positivity assumption 
has recently been discussed at length in the literature [a-llOfl. Such work is not only important for quantum error 
correction and the description of noise in quantum systems but also for quantum control more generally. (See 
Ref. [l2l fl3j for recent reviews of quantum control theory and applications.) 

Studies of maps of qubits have been extensive. Some notable discussions are found in Refs. [3, EH where they 
describe a geometric picture of the spaces of states of qubits and the ranges of maps of qubits. This is done using 
the Bloch vector, or polarization vector parameterization of the density operator jl6l - ll8| and an afiine mapping of 
the polarization vector. This afiine map is essentially equivalent to the OSR, but is more geometrical, leading to 
geometric constructions which are sometimes helpful tools for visualization as well as analysis. 

The subject of noise in quantum systems governs most of the discussion of such maps in Ref. [14], and this will 
also be the case here. The discussions in Refs. [H], EH are based on the Bloch vector parameterization of the density 
operator. The generalization of the Bloch vector is known as the generalized Bloch vector, coherence vector, or 
polarization vector. Since the dimension of this vector grows with the dimension of the Hilbert space as d 2 — 1, where 
d is the dimension of the Hilbert space, the simple three dimensional picture (for d = 2) becomes more difficult to use 
and one loses the ability to visualize the entire space. However, the polarization vector parameterization has several 
appealing properties. (1) The components of the polarization vector are measurable quantities; they are proportional 
to the expectation values of Hermitian operators. (2) The trace condition (the trace of the density operator must be 
one) and the hermiticity of the density operator become apparent in this picture. (3) The positivity of the density 
operator can be expressed in terms of the components of the density operator and the magnitude of the polarization 
vector is directly related to its purity ([l9|, H(|). (4) The Casimir invariants, quantities which are invariant under all 
unitary operators, are easily written down in terms of the polarization vector (l9| . (5) A tensor product basis may 
be used so that subsystems of interest are manifest fl7l . Il9j |. Clearly these are not entirely independent properties. 
However, they are useful properties. For these reasons, this picture, which generalizes the Bloch equations, is becoming 
more widely used in both theory and experiment. This begs the question, which we attempt to answer, how can we 
fully utilize this picture? 

In this paper we discuss the generalization of afiine maps of the polarization vector to systems with d dimensions. 
We provide explicit expressions for the affine map of the polarization vector in terms of the OSR elements. (The 
representation of a map of the polarization vector generally contains linear and translational terms.) Our explicit 
calculation provides a direct link between the OSR, pioneered by Sudarshan, Mathews and Rau [3] and the afiine map 
picture as is done in Refs. [HI, EH for qubits. We discuss the possibility of extending the description of the positivity 
domain for qubits fl5| to <i-dimensional systems and find a continuous symmetry breaking description which describes 
the inability of one to apply a simple singular value decomposition (SVD) to all density operators. This symmetry 
breaking provides a picture which compliments the work of Kimura and Kossakowski describing the positivity domains 
for density operators of d-state systems [2l[ but in a very different way. 

Specifically, in SectionHH a brief review of the derivation of the general operator-sum decomposition (OSR) from the 
dynamical map is given, along with a brief review of the polarization vector parameterization of the density operator. 
In Section [IIII the form of the affine map is derived from the OSR and some properties of the affine map are given. In 
Section ITVl physical examples of noise are given to show how the polarization vector changes for some important types 
of noise. We then examine a particularly important basis for the operators comprising the OSR in Section [V] Having 
established the relation between the OSR and affine map, we show how one can find the affine map directly from 
the corresponding dynamical map in Section I VI I Section IVIII contains a review of the singular value decomposition 
of the affine map for qubits and our attempted generalization. Finally, concluding remarks and future directions are 
provided in Section IVIII1 
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II. BACKGROUND 

In this section we provide background for what follows. This primarily follows the results of Ref. Q. 

A. Dynamical Maps and the SMR Decomposition 

As did Sudarshan, Mathews, and Rau, let us consider a general mapping from one Hermitian matrix to another of 
the form 

P > = Ap, (1) 

or more explicitly 

Pr's' — A r ' s ' rs p rs . (2) 

Since the density matrix is required to be Hermitian p — p\ positive semi-dcfinitc, p > 0, and have trace one Trp = 1 
the mapping A must have the following properties in order for it to map density operators to density operators: 

Asr,s'r' (A rsr ' s '^ , (3) 

which ensures hermiticity, 

x r x s A 3rs i r /y s i y r i ^ 0, A rr s ' r / — <5 s / r /, (4) 

which ensure positivity and the trace condition respectively. SMR then introduced a matrix _B, related to A by 
relabeling, such that 

B rr / ss ' = A sr s ' r '. (5) 

This matrix has the following properties: 

Brr',ss' — (B ss ' rr '^ 3 (6) 

corresponding to Eq. ([3]), and 

z rr r B rr ' ss i z ss ' ^ 0, B rr i rs / — S r ' s ' : (7) 

which corresponds to Eqs. @. In this article we will allow the map to be more general initially and only require the 
first of these conditions corresponding to the hermiticity condition, Eq. ([3]) or, equivalently, Eq. ©. Since B satisifies 
the hermiticity condition, it can be considered a Hermitian matrix and as such it is diagonalizable using a spectral 
decomposition, or eigenvector decomposition. Letting r/k be the eigenvalues of B the spectral decomposition can be 
written as 

B rr ,, ss ,=J2vk£k^- (8) 

k 

This is also sometimes written with indices supressed as follows: 

B = J2vkC k Cl (9) 

k 

where (Ck)rr' — Crr' are eigenvectors of B and rjk its eigenvalues. This is shown in detail in Q. The Cfc, as well as 
the density operator itself, may be written as a matrix or a vector. 

If the rjk are all positive, they may be absorbed into the C's to arrive at the familiar form of the operator-sum 
decomposition: 



>' = E^4> ( 10 ) 



where A k = ^/f^C k [Hj]. 
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B. Polarization Vector Representation of the Density Operator 



Before we present the polarization representation of density operators, it is important to provide our conventions. 
These are contained in Ref. [l9| which follow Refs. [H|,[23j]. However, they are not completely standard; see for example 
[2 G3, H3] for other conventions. A density operator on a d-dimensional Hilbert space Hd will be represented 
using a set of traceless, Hermitan matrices {A^}, i = l,2,...,eP — 1 with the normalization condition Tr(AiAj) = 
2<5.jj, commutation [Aj,Aj] = lifijk^k (fijk are the totally antisymmetric structure constants), and anticommutation 
relations {Xi,Xj} = gMy + 2dijkXk (the dijk are the totally symmetric d-tensor components), where the sum over 
repeated indices is to be understood unless otherwise stated. (In some cases the sum is displayed explicitly for 
emphasis.) These three relations may be summarized using the following equation 



2 

XiXj — ^1<% + dijkXk + ifijkXk- (11) 



The density operator can now be written as 



p=-lt + bn-\), (12) 



where b — -J (d(d — l)/2). The "dot" product is a sum over repeated indices, 

d 2 -i 

a ■ b = afii = a,ib t . (13) 
i=i 

Any complete set of d 2 — 1 mutually trace-orthogonal, Hermitian matrices can serve as a basis and can be chosen to 
satisfy the conditions given here. The components of the polarization vector n are proportional to the expectation 
values of the set of Hermitian observables A: 

ni = ^Tr(A,p). (14) 

Therefore these are directly measurable quantities which can be used to completely specify any state, pure or mixed. 
Pure states have the properties that 

n-n = l, and n-kn — n, (15) 

where the "star" product is defined by 



,~ 1 did -I) , 

{a * b)k = d^2\ 2 ( m > 

For later use, a "cross" product between two coherence vectors can also be defined by 

(a x b) k = fijkdibj. (17) 
The set of mixed states can be specified in terms of a set of positivity conditions [l9l . [20| . 



III. AFFINE MAPS FROM THE OSR 



In this section we obtain the affine map of the polarization vector in terms of the components of the OSR of the 
dynamical map. However, we begin with a general case, as is done in Ref. [15j . and restrict to particular classes of 
maps which are often physically relevant. Furthermore, we do not restrict to completely positive maps. The affine 
map provides a final connection between these three different forms of the map: the dynamical map A (or equivalently 
B), the OSR, and the affine map. 



A. Explicit form for the Affine Map 

Let the operator-sum decomposition for a map from one density operator to another be given by 

p' = $(p) =J2vkC kP Cl (18) 



Following Ref. 15] which treats the qubit case, the C k can be represented by (complex) linear combinations of the 
{Xi}, as 



C k = v Qk l + v k ■ A, 

where v ok , v ik G C. 

The following identity will be used repeatedly in our derivations 

(at + u ■ X)(bl + w • A) = ^ab + — u ■ w^j 1 + ^aw + bu + iu x w + —u-kw j • A. 

where the definitions Eqs. (| 10|) and (|17p were used along with the following definition, 



d(d-l) 1 



d-2 



(19) 



(20) 



(21) 



Two special classes of maps are particularly important, unital maps and trace preserving maps. They are defined 

by 

1. $ is called unital if $(1) = 1, 

2. $ is called trace preserving if V p, Tr(<f>(p)) = Tr(p). 

If the map is unital, they by explicit calculation using Eq. (1181) . 

2 



E 



Vk 



\v ok \ +-;v k -v* k 



= 1, 



and 



E ? ? fc ^ofeWfe + wS/c^fc + wfe x ^ + i« fc * ^ = 0. 



Similarly, if the map is trace preserving, then an explicit calculation shows that 

2. 



^2vk[\vok\ z + ^v k -v* k 



1, 



(22) 



(23) 



(24) 



and 



»7fc ^UofcWfe + VQ k v k +iv* k xv k + ^-v* k * v k ^j = 0. 



(25) 



Note that if the map is either unital or trace-preserving, then the condition given in Eq. (|22|) (equivalently Eq. (|24|) ) 
holds. 

To provide the explicit form of the map, we first use Eqs. (fT8|) , (fT9|) . and (|20|) . so that 

(26) 



*(p) = E = E Vk(vokt+ v k -\)-(l + bn- \)(v* ok l + v* k ■ A), 



can be rewritten and the identity and A parts treated separately. 

First, consider the identity part. The coefficient of the identity can be written in the following form 



E% 



\V0k\ 



+ ^ ( v o k Vk + vokV* k +iv* k x v k + h% -k v k ] • /7 ) . 



(27) 
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If the map is unital, using Eqs. (|22j) and (|23|) . the coefficient of the identity reduces to 

1 Mb 



d <P 



J2vk{v* k x v k ) 



(28) 



If the map is also trace-perserving, subtracting Eq. (|23|) from Eq. (|25|) implies that the coefficient of the identity is 
l/d. 

Now let us consider the non-identity part of the map. Denoting the result of the map by p' , we have 



=p' = -(l + bn'-X). 



Viewing the map as an affine map of the coherence vector n, 

«H>n' = Tn + £, 

so the components of ft' are given by 



Tlq TpqTlp -f- tq. 



Thus, after some rearrangement, the map can be written as 



-y 



2b. 



Vk 



1 



v ok v k + v ok v k + iv k x v k + ~v k * v k 



■A 



(29) 

(30) 
(31) 

(32) 



% fc| n+ — {v k ■ n)v* k + ib{vg k v k x n + v ok n x + -(v ok n-kvl + Vg k n*v k ) 



+—{{vk *ri) x vt + (v k x n) *vt\ - 6(i7fc x n) x vt + \(v k -kn)*vt 
c c 



■A, 



where T = 52 fc 77 fe T fe and i = £ fc The fc th term <J3U), is 



1 



v ok v k + v ak v k + iv k x v k + -v k ★ i> fe 



(33) 



(34) 



Note that t is zero for a unital map. 

The k th term represents the action of the real (d 2 — 1) x (cP — 1) matrix T. (The proof of the reality of T is 
given below.) The components of the linear part of the map, which we refer to as the T matrix, may then be written 
as 



Vk 



\v k\ - jV*k v rk ) S PQ + -^( v pkV* qk + V*kV q k) + ifrpq(VQ k V rk ~ V 0k V* k ) 



+ d rpq (vo k v rk + v Qk v* k ) - id rpq f srt v* tk v sk 

{dqsrdrtp ^qtv^spr &qpr3>str)^tkP sk 



(35) 



It is worth emphasizing that the T matrix and the vector t specify an affine map of the polarization vector. We also 
note that the map is linear when t = 0. This is equivalent to the dynamical map, A (or B) specified in Section HIl and 
also the operator-sum decomposition. Thus we have provided a mapping of a vector to a vector which corresponds 
to a general mapping of a density operator to a density operator. The next question which we will answer in part is 
the following. What are the properties of this map? We are able to provide partial answers to this broad question for 
both the general map and for some particular cases of interest, which we do next. 



B. Properties of the Associated Linear Map 



In this subsection we show that the linear term in the affine map, i.e., the T matrix, is real. Although this must 
be true, it is shown explicitly since it is not obvious from the expression Eq. (I35[) . In addition, we provide general 
conditions for the matrix to be symmetric. 
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1. The T '-matrix is real 

It is not clear that the third to the last term in the curly brackets ({}), ([33]) . which mixes star and cross products, 
is real. If we can show that this is real, then we will have shown that T is real since the other terms are obviously 
real. Let us consider whether 

+ — X! Vk [((4 *n) X v%) ■ A + ((v k x n) * v* k ) ■ A] (36) 
c k 

is equal to its complex conjugate. 

To provide a sufficient condition, let us first consider a term with a given k. From the identity, Eq. (|A3|) : 

dijm fmln ~\~ djlmfmin ~\~ dlimfmjn — (^^) 

we can form the following identities with any vectors a, b, c, A. (Here we will use A as the basis for the Lie algebra, 
but in these identities, it could be any vector.) First, we contract the LHS of Eq. (I37|) with aibjCi\ n and find: 

[(a-kb) x c] ■ A + [(b-kc) x a] ■ A + [(c-ka) x b] ■ A = 0. (38) 
For the next identity, we contract Eq. Q37[) with ajbic n \i to obtain: 

[(6 x c) * a] • A + A • [c x (a * 6)] + [(a x c) * 6] • A = 0. (39) 

For our purposes, it is relevant to note the following symmetries in the indices. The identity is invariant under the 
interchange of the following pairs of indices (I, i), so that these are the only two different types of identities 

when we distinguish A as a set of basis elements. 

We will use this second identity Eq. (|3T)1) to rewrite the first of the two terms in Eq. (|3^|) . First, note that 

[tffc x (v k -kn)] ■ A + [(n x u£) • A + \(v k x u£) *n] • A = 0. (40) 

Now the first term in this equation is the same (up to sign) as the first of the two terms in Eq. (|36l) . This allows us 
to rewrite Eq. (|36|) as 



ib 
c 



|[(n x vl) * v k ] ■ X + [(vk x vl) -kn] ■ X+ [(vu x n) *v* k ] ■ X| . (41) 



We now want to show that the complex conjugate of the coefficient of A is real by showing it is equal to itself 
(z* = z z is real). Let us take the complex conjugate of the coefficient recalling that n is real: 

^ { [(n x v k ) * v* k ] ■ A + [(v* k x v k ) * n] ■ A + [(v* k x n) * v k ] ■ X} . (42) 

The first term here is the negative of the third term above, the second is the negative of the second above and the 
third is the negative of the first. Therefore the two are equal due to the overall minus sign in this latter expression 
and so the coefficient is real and this part of the map is real. To prove that this sufficient condition is also necessary, 
we use the fact that there exists a minimal decompostion of the map such that each T k is independent and can always 
be put in minimal form [24]. Then, since the r\ k a re real, this condition is also necessary when the map is in minimal 
form. □ 

2. Conditions for T to be a symmetric matrix 
The importance of this stems from the fact that a symmetric matrix can be decomposed in a polar decomposition 

M 

S = ODO, (43) 

where S is a real symmetric matrix and O is orthogonal. We will show this to be true for a given k and the 
generalization follows from the sum J2k ^k- 
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The only antisymmetric term in the expression for T is 

T pq = Vkifr P q(vo k v rk - v ok v* k ) (44) 
k 

Therefore, for T to be symmetric, we require this term to vanish. 

Again, let us first consider the necessary condition that each term given by a fixed k vanish independently. Since 
the matrices (faff)-/ form a representation of the Lie algebra of SU(d), the only way for this to to happen is if all 
vector components, indexed by 7, vanish independently. This implies that we must have 

VokVyk = v kv* k . (45) 

There are several ways in which this can happen. Let us consider some examples. Any of the following is sufficient to 
ensure that T is symmetric: 

1. vok and v lk are all real. (Thus all C k are Hermitian.) 

2. For a given k either VQ k = 0, or v lk = 0, V7. (Note v lk — for all 7 is rather trivial because then C k oc 1.) 

3. v lk — vg k , Vfc, V7,<5 and vo k = v jk . Note that the vectors can be different for each k, but all have the same 
components for each k. They are also NOT required to be real. 

IV. EXAMPLE CHANNELS 

Here we consider some examples of maps which are of special interest. These are far from exhaustive and we will 
consider other examples in future applications. 

A. Ck hermitian 

When the C k are hermitian, the T matrix is symmetric as shown in the previous section. Here we provide explicit 
expressions. 

Consider a map with C k = C\ so that the C k are Hermitian, or sometimes called self-adjoint. Then C k can be 
expressed as 

C k = (vok^ + VikXi), (46) 

with v ok ,v ik £ K. 

Since the components of T, Eq. ()35|) . are real, the matrix components of the affine map for Hermitian C k , are 

2 2 

(T k ) pq = |«ofc| S pq + -(v k ) p (v k ) q + 2v ok d rpq (v k ) r 

-ftprfrsq(v k )t(v k ) s + dtprd r sq(v k )t(v k ) s , (47) 

where the third and fifth terms in Eq. (|35[) vanish since v o k , v k are real. This is symmetric (in p and q) since the first 
two terms clearly are symmetric and the latter two can be shown to be symmetric by renaming the dummy indices s 
and t, and cyclicly permuting the indices on c and d, viz. 

ftprfrsq(v k ) t {v k ) s = f spr frtq(v k ) s (v k )t 
= frspftqr(v k ) s (v k )t 

= ftqrfrsp(v k )t(v k )s- (48) 

Therefore if the C k are Hermitian, then the map T pq is real and symmetric in p and q. 
The affine part of the transformation reduces to 

t k = \{2v Qk v k + -v k *v k ). (49) 
c 
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B. C k Unitary 

This case will be explored further in the next subsection where numerous reasons are given for the special consid- 
eration. Here we provide an explicit expression for the T matrix and note that t is zero for this case. One case of 
interest is when the map is given by $(p) = J2 k r lkU k pU k where J2 k ^ k = an< ^ ^ fe > 0- 

Consider a map with C k C'l = 1 = C'lC'k. Then C k can be expressed as 



with the following constraints 



and 



This implies separately that 



and 



This implies that T k reduces to 



C k = {v ok t + v lk \i), (50) 



\v 0k \ 2 + ^v k -v* k ) I. (.VI i 



vokV* k + v* ok v k + iv k x vt + - v k * v* k 

c 

VOkVl + VQ k V k + iv* k X V k + ^V k -kVk) ( 52 ) 



v ok v* k + v* ok v k + ~v k *v* k ) ( 58 i 



= v k xv* k . (54) 



2 

{T k ) pq = \v 0k \ 2 S pq + ^(Vk)p(^)q + i[Vo k f rpq (Vk)r ~ V kfrpq(Vk)r] 

(v k ) r j + i[d S prfrtq(vk)s(Vk)t + d r sqftpr{v k )t{v* k )s\ 
-ftprfrsq(vk)t(Vk)s + d tp rd rsq (v k )t{v* k ) s . (55) 



C. Qutrits 



In this section we provide explicit example channels which are relevant for physical problems involving three-state 
systems. See also Ref. [25[ and references therein. 



1. Qutrit Depolarizing Channel 



In general the effect of the depolarizing channel on a density operator is to cause a uniform shrinking of the 
polarization vector ft. This is written as 

p = -(t + bn- A) -> p' = -(t + bpft-X), (56) 
d d 

where p is the shrinking factor. (See for example [26| and references therein.) 

The effect of the depolarizing channel on the coherence vector is such that it uniformly shrinks each component of 
ft by a common factor. It can be readily verified that if we let 

Co = Vl~x 1, (57) 

and 



C k = ^3x/16X k , fc=l,2,...,8 



(58) 
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we then have J2k=o C kC' k = 1, t = 0, (T ) pq = (1 - ir)<W and ( T fc) P9 = {x/8)6 kp 5 kq + (3x/16)(d kpr d rkq - fkprfrkq) 
for fc = 1, 2, . . . , 8. The coherence vector is thus transformed according to 

n^Tn = ^T k n = (1 - 9x/8)n. (59) 
Thus the shrinking factor is p = 1 — 9x/8. 



2. Qutrit Phase Damping Channel 



Analogous to the qubit phase damping channel, the effect of the phase damping channel on a three state system is 
such that it leaves the diagonal components of n unchanged while uniformily shrinking the off-diagonal components. 
When the operators C k are chosen to be 

(60) 

(61) 
and 

C 3 = V3x/20(1 + A 8 ), C 4 = V3x/20(1- A 8 ), (62) 

it can be shown that YX=o C kC\ = 1, t = 0, (T ) pq = (1 - x)S pq , {Ti) pq = fi£, (T 2 ) pq = M3 , { T s)pq = Ms > and 
(T 4 ) pg = yUg", where 



Co = 1, 
d = y/3x/20(l + X 3 ), C 2 = ^3^/20(1- A 3 ), 



M; 



± _ 



3a; 
20 



r^pl^ql dl p q flprflqr + ^Ipv^lqr 



(63) 



When r] k = 1 for Vfc the matrix T shrinks the off-diagonal components of the coherence vector by a factor of 1 — 3x/5, 
i.e., 



rii, for i = 3, ! 

(1 — 3x/5)nj, otherwise 



(64) 



In the last expression we have taken the usual convention of labeling the diagonal elements of the basis as 3 and 8. 



3. Off-Diagonal Channel 

Let us now define 

Co = 1, (65) 

and 

C k = ^/x/l\ k , for k = 1,2,4,5,6,7. (66) 

These values of k correspond to the off-diagonal generators of SU(3). It can be verified that this particular OSR 
transforms the coherence vector according to 

^ f (l-3x/2H, for i = 3,8 m 
1 \ (1 — x)rii, otherwise v ' 

when r\ k = 1 for Vfc. In this case, the diagonal components of n shrink faster with x than the off-diagonal components. 
If instead r\ k = (-l) fe , the mapping given by T = J2 k (~^-) k Tk (again t = Y^ k (~^-) k tk = 0) affects the components of 
n in such a way as to shrink them according to 

( (1 - 3x/2)m, for i = 1,4,6 
n>i 1 ^ \ (l-x/2)m, fori = 2, 5, 7 (68) 
I (1 — x)rii, for i = 3, 8 

The components of n associated with the real off-diagonal elements (Aj, i = 1,4,6) shrink with x faster than 
those associated with the imaginary off-diagonal elements (A;, i = 2, 5, 7). The diagonal components are reduced in 
magnitude at an intermediate rate. The values assigned to the i] k have a significant effect on the evolved state, as 
clearly illustrated by this example. 
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4- Trit Flip Channel 

For qubits, the bit flip channel acts in such a way as to flip the basis state |1) to |2) and vice versa. For qutrits 
there are three basis states which can be labeled as |1), |2), and |3). To describe the effect of flipping one of these 
basis states to another in terms of an affine mapping of the polarization vector let us first define the following pure 
state density matrices 



Pi 





P 2 = ( 10 





P3=\ ooo 
1 



(69) 



Now the affect of flipping from the state |1) to |2) is equivalent to the density operator p\ evolving to p 2 . 
be acheived by the transformation 2 )P^\\ 2 ) ~ <° 2 > wnere 



This can 



2) 




(70) 



Notice that C^ 2 ) is an element of the Gell-Mann basis for SU(3). The affine map associated with this transformation 
has a linear part as well as a nonzero translational part, specifically 



T (l,2) 
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-1 


























-1 

















































































































\o 




















1/3 J 



and t( 1;2 ) 



/ \ 









V 1/3/ 



(71) 



Note that since 2 )P^h 2) = P^i * ne amne ma P associated with the flip |2) — > 1 1) is the same as above. Similarly 



one can show that the transformation 3 )PiC( 1 3 ) = P3, with 



C(l,3) 




(72) 



is equivalent to the affine mapping 



T, 



(1,3) 



/0 













Vo 











-l/Vs 









10 
0-100 














-1/V3 





-2/3 / 



, and i (li3 ) 






1/2^/3 





V -1/6 J 



(73) 



Since C, x 3^3^! 3 -. = pi, the mapping associated with and ^1,3- 



will yield the transition |3) — > |1). For the 



flips 1 2) <H> 1 3), the associated OSR decomposition has the term C( 2 ,3) with 



(2,3) 




(74) 
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This leads to an affine mapping of the form 



Ti 



(2,3) 



/0 







Vo 





1 



0-1 



l/y/3 



\ 



i/V5 






-2/3/ 



and t 



(2,3) 






-1/2V3 





V -1/6 J 



(75) 



V. A PREFERRED BASIS: C k UNITARY AND ORTHOGONAL 

It is often the case that a preferred basis for the Ck is chosen. One important case, is when the Ck are both orthogonal 
and unitary. This choice corresponds to a nice error basis [27| and applies to the so-called generalized Pauli matrices 
which are formed from Ck given by U mn = X m Z n , m, n = 0, 1, 2 where X\j) = \j + lmodd), Z\j) = u^\j), w — g 2 W3 
and {|i)}f=o is an orthonormal basis. We also note that it can be applied to decoupling controls which form an 
discrete subgroup of the unitary group, but in that case the Hamiltonian is modified rather than the density operator 

& 

Since the case where the Ck are unitary and orthogonal is a particularly important case, we show that when the 
Ck are orthogonal, unitary, and belong to a discrete subgroup of the unitary group, the trace of pairs of unequal Tk 
must be —1 between elements of the set. 

Let {Uk} be such a set which also forms a subset of a fundamental irreducible representation of a discrete subgroup 
of the unitary group, i.e., UjjJk = % and the set of Uk can be extended to form a discrete subgroup of the unitary 

group, and Tr(UiUj) = for i ^ j. Then, since UjUi is a group element, the trace corresponds to the character of 
the product. The character, for an irreducible representation, is the same for each element in a class and, in this 
case, the class has character equal to zero. However, it can be different for different representations. Note that the 
corresponding T matrix is an orthogonal matrix since it belongs to the adjoint representation of the group which is a 
real (d — l)-dimensional representation. This follows from the equation 

UkpUl = [/ fc i (l + brl • A) u\ = J f 1 + bJ^n i (R k ) ij X^j , (76) 

where Uk € U(d) and Rk is in the adjoint representation of the group. In this case Tk corresponds to Rk, and Uk is 
the same group element as Rk but in a different representation. 

Since Tr(UjUj) = 0, and Ui (§5 Uj is in the tensor product of the representation of U(d) ® U(d), we can write 
Tr([// (g> UjUk ® U{) — and this is the product of the two representations. If we take the complex conjugate of the 
first of these two representations (to obtain the conjugate representation) we find that the product, which is equal to 
the direct sum of the real (d 2 — l)-dimensional representation and the trivial one, must be zero. Since the trace of 
the sum is the sum of the traces, and the trivial representation has character (trace) one, the other must have trace 
— 1. Therefore, if the Ck are unitary, orthogonal, and the set of these can be extended to form a group, then the 
corresponding Tk have the traces of the products all equal to —1, viz., Tr(CiCj) = <f> Tr(T^Tj) = —1, V i ^ j. 

As mentioned above, this result is important for error prevention and control methods. It is also generally important 
for translating between the OSR and the affine map when this basis is chosen. 

We now briefly discuss how one would obtain the affine map directly from the dynamical map, showing that the 
OSR is in fact not a necessary intermediate step. 



VI. AFFINE MAP FROM THE DYNAMICAL MAP 



It is fairly straight-forward to show that the affine map can be obtained directly from the dynamical map. In this 
case, a convenient basis is the set of matrices which form a basis for GL(n, C). Such a basis is defined by the set of 
matrices with a 1 in the ith column and jth row and zeros everywhere else, denoted {Eij}. Given that from Eq. @,we 
see that the two bases may be converted from one to another using, for example, the correspodence between the raising 
and lower operators of the group and the hermitian basis. To be specific, one would expand p = s p rs E rs , and 
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use the corresponding hermitian basis A rs — E rs + E sr and A' rs = i(E rs — E sr ) as well as linear combinations of the 
identity and Cartan algebra elements to expand the diagonal matrices. (See for example (29| and references therein.) 

Thus we have given the following correspondences: dynamical maps — > OSR — > affine map and dynamical map — > 
afhnc map. Note, however, that the OSR is not unique - not for completely positive matrices 14] nor for those that 
are not completely positive 24]. Therefore, the set of Tk matrices is also not unique and thus converting from one 
picture to another is not a one-to-one transformation. The dynamical map is essentially unique as is the T matrix as 
well as the minimal decomposition of the map (36j . 



VII. SINGULAR VALUE DECOMPOSITION 



In the analysis of maps of qubit density operators, sufficient conditions for the positivity of maps follows from 
positivity conditions described by the Bloch-sphere (la. [30j. The maps, represented as an affine map of the vector, 
can be decomposed in terms of a singular value decomposition (SVD) which greatly simplifies the analysis Here 
we discuss a decomposition which is similar to the singular value decomposition of the affine map for qubits, but 
applies to a <i-dimensional system when d > 2. This turns out to be much more complicated for d-state systems than 
it is for a two-state system. The reason, as explained in detail below, is that there are restrictions on the "rotation" 
part of the map (that part which does not change the magnitude of ft). The restrictions are determined by the 
"shrinking" part of the map (that part which can reduce the magnitude of one or more components of ft) and differ 
for different shrinking matrices. Thus we are seeking an answer to the following question. To what extent can we 
generalize the results of qubits to qudits with d > 2? Our results provide conditions under which a SVD may be 
performed as well as restrictions to its use. This provides some insight into the geometry of the space. 



A. SVD For Affine Qubit Maps 

Before presenting our results for qudits, we will review the basic idea for qubits as used in Ref. [lEj]. Let some initial 
density matrix 

p=\{l + n-5), (77) 



be acted on by 



$(p) = l(l + n'-a). (78) 



This map can be described by an affine map of the polarization vector, 

ft' = Tn + t. (79) 

In what follows we will let t — and discuss T. The translation may be treated separately. 

For two-state systems, maps can be diagonalized by means of a SVD. This is because any real N x N matrix M 
can be written as 

M = 1 D0 2) (80) 

where 0, € SO(N) and D is a diagonal matrix with the entries being called the singular values. For qubit maps 
TV = d 2 — 1 is 3 and an action of £7 € £7(2) on the density operator is equivalent to an SO(3) matrix acting on the 
polarization vector ft. To see this, we write 



(81) 



where R G 50(3). Thus acting before and after a map enables one to choose a preferred basis for the map. This 
basis can be chosen so that T is diagonal since choosing R and R' appropriately enables us to put T into a diagonal 
form using D = RTR' . This set of diagonalized maps can be considered a double coset space of the elements of the 
set of maps. 

For future reference it is relevant to note that for a positive map, the matrix D cannot increase the magnitude of 
ft or any component of fx. (Again, we are taking t — 0.) In general, for the SVDs we will consider here, we will refer 
to the D part (diagonal part) as the "shrinking factor" since, for a positive map, it can decrease, but not increase, 
the magnitude of the polarization vector. The parts O will be called "rotations" although the maps can clearly be 
physically achieved in many different ways. 
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B. SVD For Qutrits 

Here we will discuss the case of qutrits. The generalization to qudits follows. We know that for a mapping from 
pure states to pure states, the map is a unitary transformation. That is, any map from a pure state to a pure state can 
be written as a unitary transformation. The unitary transformation acts on the density operator to produce an affine 
map which has the form T eAd(SU(3))cSO(8). That is, the T matrix is an element of the adjoint representation of 
SU(3) which is a subset of SO(8). It clearly does not span the space of SO(8) in the way that SU(2) covers SO(3) 
since the dimension of SU(3) is 3 2 — 1 = 8 and the dimension of SO(8) is 8 ■ 7/2 = 28. 

For an even simpler case, consider a map which acts as a depolarizing channel together with a "rotation." In this 
case there is a uniform shrinking factor such that each component of the polarization vector is reduced by the same 
factor p < 1. (See for example [26| and references therein.) For qutrits, when p < 1/2, the map has a SVD with 

T = 1 D0 2 , (82) 

where Oi € SO(8). Any matrix in SO(8) is allowed since any direction for the polarization vector leads to a positive 
density matrix out put as long as the shrinking factor is p < 1/2 and uniform. This follows from the positivity 
conditions in Refs. (Il|2(j. 

For there to exist a SVD for maps which have 1/2 < p < 1 the SVD can not either be the SVD for p < 1/2 nor 
can it be the SVD for pure states. The SVD for these intermediate values of p must have some elements Oi € SO (8) 
that can not be arbitrary elements of SO(8), but come from a restricted set which depend on the shrinking factor. 
However, they are also not restriced only to the set of matrices in SU(3)c SO(8), they are from a larger set. The 
question is, what is that set and how do we find the values? 

Whereas the general problem is quite complicated, we can make some general statements about the allowable 
rotations. Let us consider the algebra of SU(3)c SO(8). For pure states, the set of allowable rotations are only 
those which are in SU(3) and at the other extreme, all rotations in SO (8) are allowable when the magnitude of the 
polarization vector |n| < 1/2 as discussed above. Thus the SO(8) symmetry is much larger than the SU(3) symmetry 
and it could be said that the symmetry is broken continuously by the increasing magnitude of the polarization vector. 
Such considerations are important for a variety of reasons. Most notably, the "good quantum numbers" are labels 
for states which are carriers of an irreducible representation of a group. (See for example [3lj.) Thus symmetries 
determine simultaneously measurable quantities for a quantum system and thus enable us to specify a state as well as 
we can. When a symmetry is broken, the symmetry for the system is reduced and the set of good quantum numbers 
changes, indicating a change in allowed states of the system. In this description of the allowable states, one possible 
symmetry breaking construction which could be used for broken symmetries is the following subgroup chain: 

50(8) D SO{6) x S0(2) = SU(4)/Z 2 x 17(1) D SU(3), 

Note, however, that the restricted symmetry arises from positivity conditions on the density operator. This positivity 
condition is associated with purity/mixedness of the system. Thus we have a connection between the positivity of 
the system and its entropy. This relation is only direct for two-state systems, however, since there are two positivity 
constraints for three-level systems and more for higher-dimensional systems. 

The explicit parameterization of the rotation matrices can be obtained from the exponentiation of the algebraic 
elements given in the appendix. There we have expressed the algebraic elements of S0(8) in terms of the algebraic 
elements of SU(3). This embedding of the SU(3) rotations into the set of S0(8) rotations is important since it specifies 
the unitary part of the transformation. We will not provide the explicit parameterizations in each of the regimes, but 
do discuss an example. 

The general case is straight-forward to infer. There is a ball at the center of the space of density operators which 
has a spherical symmetry in which all density operators are positive semi-definite. At the other extreme are the pure 
states which can only be acted upon by a rotation part consisting of unitary transformations. The intermediate points 
have "rotation" matrices which are restricted and depend on the magnitude and direction of the polarization vector. 

As an illustrative example we consider the case where a qutrit is initially in the pure state |2) (let |1),|2), and |3) 
represent the basis states of a three-state system). The coherence vector associated with the density matrix |2)(2| is 
given by n = (0, 0, — V3/2, 0, 0, 0, 0, 1/2) T , where the third and eighth components are associated with the diagonal 
generators of SU (3), i.e., |1)(1| — 12)(2| and l/v / 3(|l)(l| + 12)(2| — 2 |3)(3|) respectively. If the components of n experience 
a rotation due to, say, the orthogonal matrix f 2 s (see Appendix B), and shrink in magnitude by the factor (1 —p), 
the evolved state nf = (1 — p) cxp {—i6f2%)n may or may not correspond to a positive-semidefinite matrix. Since the 
initial state was pure, the necessary condition for positivity is S3 = 1 — 3ri' ■ n' + 2(n' * n') ■ ft' > 0, see Ref. (l9l. [20|. 

In Fig. [T] we plot the region of positivity as a function of the parameters 6 and p. The figure shows that while S3 — 
initially (p = 6 = 0), this quantity quickly becomes negative as the state begins to rotate without shrinking. The 
figure also shows, although not very clearly, how the initial state becomes positive again as 9 increases to multiples of 
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FIG. 1: Plot showing the region of positivity when the initial pure state \2) is allowed to experience a uniform shrinking (1 — p) 
and rotation of angle 9 by the matrix /28 (see Appendix B). The quantity 1 — 3n' ■ ft' + 2 (ft' *ft') ■ ft' is graphed vertically since 
for this initial pure state the condition 1 — 3ft' ■ ft' + 2(ft' * ft') ■ ft' > is necessary and sufficient for positivity. 

27r/3 while no shrinking has occured. When the components of the initial coherence vector uniformily shrink by an 
amount (1 — p) > 1/2 we see that any rotation about the matrix /2s leads to a positive density operator. 

VIII. CONCLUDING REMARKS 

We have obtained an expression for an arbitrary affine mapping of the polarization vector associated with a d- 
dimensional density operator in terms of the components of the operator-sum representation of the dynamical map. 
This provides a geometric picture, albeit somewhat abstract, via the polarization vector components which are mea- 
surable quantities. We also described a direct method for expressing the affine map in terms of the dynamical map. 
Some example channels were provided in order to highlight the connection between particular terms in our expression 
and specific operators appearing in the OSR. For the important case that the OSR components are unitary and 
orthogonal, we have shown that the corresponding affine components are real and have trace -1 between pairs of 
unequal terms. This is particularly useful for the generalization of the Bloch equations since this basis is used for 
error prevention methods. 

To generalize the methods used for qubit maps, we discussed the possible generalization of the singular value 
decomposition of qubit affine maps to higher dimensional systems. We have found that a symmetry-breaking occurs 
due to positivity constraints. Thus we have given a relation between the physical symmetry, the purity and the 
entropy of the physical system. In this particular case we have provided the example of the depolarizing channel, but 
the generalization is also discussed. 

While we have provided some particular expressions for physically motivated quantum channels (maps) , applications 
of this work range from information theory to practical experiments as is the case with the OSR. One example of the 
utility is given in Ref. [32j . where the robust simulation of quantum systems with quantum systems was considered. 
In addition, the preferred basis given here provides a connection between quantum error correcting codes and direct 
application thereof. This is due, in part, to the ability to more directly determine the input and output, and therefore 
the map itself, for quantum systems undergoing noisy evolution. In future work we will use this to try to provide 
more practical methods for determining the best combination of error prevention methods. (See [33j j for a review.) 

More generally, we anticipate applications of this could be more wide spread in quantum control. For example, 
classical control systems utilizing affine maps may be found in Judjevic's book, "Geometric Control Theory," chapter 
4. This includes Linear systems reachability, and constraints. We also hope that our work will be useful in ways not 
yet foreseen by us. 

Appendix A: Identities 

In this appendix we provide identities that were used in the derivations of the expressions above. 

The commutation relations, anti-commutation relations, and normalization of the matrices representing the basis 
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for the Lie algebra can be summarized by the following equation: 

2 

AiAj = —Sij + ifijk^k + dijkXk, (Al) 

where here, and throughout this appendix, a sum over repeated indices is understood. 
As with any Lie algebra we have the Jacobi identity: 

filmfjkl + fjlmfkil + fklmfijl = 0. (A2) 

There is also a Jacobi-like identity, 

filmdjkl + fjlmdkil + fklmdijl = 0, (A3) 

which was given by Macfarlane, et al. [Hj]. 

The following identities, are also provided in (34|, 

d llk = 0, (A4) 

dijkfijk = 0, (A5) 

fijkfijk = dSu, (A6) 
d 2 - 4 

dijkdijk = — j — Su, (A7) 



and 



and finally 



fijmfklm — -j(fiikfijl ~ ^il^jk) + (dikm,dji m — dj km du m ) (A8) 



f piq f qj rfrk P = (J) hk, (A9) 

dpiqfqjrfrkp (-^10) 

dpiqdqjrfrkp = ( ^ J /ijfe, (AH) 
fd 2 - 12 \ 

dpiqdqjrdrkp — ( 1 dijk- (A12) 

The proofs of these are fairly straight-forward and are omitted. 

Appendix B: Algebra of rotation matrices 

In this section we provide the basis elements for the symmetry breaking from the SO(8) symmetry to the SU(3) 
symmetry. A basis for the algebra of SO (8) is given by one set of matrices and the SU(3) algebraic basis elements are 
expressed in terms of that basis using the structure constants which provide the algebra for the adjoint representation 
of the group. This provides the afforementioned embedding of the group SU(3) into SO(8) by the exponentiation of 
the algebra. 

The structure constants for SU(3) are fyk and we write them in matrix form as {fi)jk, i-e. the matrix /j has 
elements j, k. Consider the basis for a matrix representation of SO(8) given by the set of all 

m jk =i{\m-\m\i 



for which j < k — 2, 8. 
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For SU(3) we label the matrices 1-8, with elements j, k corresponding to the eight Gell-Mann matrices. The basis 



is 



fijk = i 
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These are clearly basis elements for the adjoint representation of the algebra which is a subset of the SO (8) algebra. 

We now make the connection between the two by writing the basis elements for SU(3) in terms of the SO(8) basis 
elements. 



h 


= rn 2 3 + i (m 47 


- m 56 ), 




h 




h 


= m 12 + i(m 45 


- m 67 ), 




h 


= - 2^ mi7 + 71226 + 17135 > + ~ m5S 


h 


= ^{mi6-m 2 7 


+ m 34 ) - 


— ^48, 


h 


1 . , V3 
= -^(-m 15 + m 24 + m 37 ) + — m 7 $ 


h 


= ^(mu + m 2 5 


- m 36 ) - 




h 


= "^( m 45 +m 67 ) 



(Bl) 
(B2) 
(B3) 
(B4) 



Let us now complete the basis by finding a complete set of 28 matrices which are orthogonal in the Hilbert-Schmidt 
sense and have fi, i = 1, ...,8 as members. Considering the Hilbert-Schmidt to be an inner product, we can immedi- 
ately write down a complete set. We will number them in no particular order. (IMPORTANT: We will also assume 
for now that they are normalized such that Tr(mjjmfci) = SikSji. This can be adjusted later if it is not true.) Let us 
consider a set of matrices which span the space TO23, TO47, m 56 when /1 is included. The set is (for fx) 



h = ^("^47 + m m ), f w 



-(m 47 - m 56 ). 
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Similarly, we will group the remaining matrices: (/2) 

fu = TO13 + (m i6 + m 57 ), /12 = (m 46 - m 57 ), 

(/ 3 and / 8 ) 

/13 = mia - (m 45 - m 67 ), 

(/*) 

/14 = (mi7 + m 2 6 + "135) + V3 m 58 , /15 

(/ 5 ) 

/17 = (mi6 - to 2 7 + m 34 ) + V3 m 4 8, /is 

(/e) 

/20 = (— TO15 + TO 24 + TO37) - V3 m 78 , /21 

(/7) 

/23 = (TOl4 + m 2 5-TO36) + V / 3m 68 , / 2 4 : 

Finally, we have the following three: 

/26 = Tils, /27 = m 2 8, /28 = "1 38 . 
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